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Abstract 

A survey  is  given  of  algorithms  for  passing  a curve  through  data  points  so  as  to  preserve 
the  shape  of  the  data. 

1 Introduction 

We  consider  the  problem  of  passing  a curve  through  a finite  sequence  of  points.  We  want 
the  curve  to  preserve  in  some  sense  the  shape  of  the  data,  i.e.  the  shape  of  the  curve 
gained  by  joining  the  data  by  straight  line  segments  (which  we  call  the  ‘piecewise  linear 
interpolant’).  We  do  not  consider  the  important  problems  of  approximating  the  data 
by  a curve,  or  of  shape-preserving  interpolation  by  a surface.  The  short  length  of  the 
paper  forces  it  to  be  selective.  So  we  concentrate  on  actual  algorithms  for  solving  the 
problem  rather  than  related  theory.  Also  we  consider  only  algorithms  where  the  curve 
is  defined  explicitly,  not  implicitly  either  as  the  zero  set  of  a function  or  as  the  limit  of 
a subdivision  process  (though  there  are,  to  our  knowledge,  extremely  few  such  implicit 
shape-preserving  schemes). 

In  Section  2,  we  consider  planar  curves  given  by  a function  y = f(x),  often  rather 
misleadingly  referred  to  as  ‘functional  interpolation’.  There  are  numerous  such  schemes, 
dating  from  1966,  with  most  of  them  prior  to  1990.  Our  treatment  is  therefore  very 
selective.  Section  3 deals  with  parametrically  defined  planar  curves,  for  which  the  schemes 
are  fewer  and  more  recent.  Finally,  in  Section  4,  we  consider  curves  in  three  dimensions, 
often  called  ‘space  curves’.  Here  the  work  is  much  more  limited,  dating  only  from  1997. 

We  note  that  in  shape-preserving  interpolation,  the  map  from  the  data  to  the  function 
describing  the  curve  must  be  non-linear.  In  what  we  call  ‘tension  methods’  the  curve 
can  be  constructed  by  a linear  scheme  for  any  choice  of  certain  ‘tension  parameters’. 
These  parameters  are  then  varied  so  as  to  ‘pull’  the  curve  towards  the  piecewise  linear 
interpolant  until  the  shape  criteria  are  satisfied.  Though  there  are  a few  variations  on 
this  theme,  there  is  generally  a clear  distinction  between  tension  methods  and  other 
schemes,  which  we  shall  term  ‘direct  methods’. 

2 Functional  interpolation 

Given  data 

(: Xi,yi)  € Ft 2,  i = 0,...,N,  x0  < xi  < ■ ■ ■ < xn,  (2.1) 
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we  consider  a function  /:  [xo,  x jy]  — > R satisfying 

f{xi)  — Vi:  i = 0,...,N.  (2.2) 

For  some  reasons,  perhaps  the  physical  situation  which  / is  intended  to  model,  we 
may  wish  the  graph  of  / to  inherit  certain  shape  properties  of  the  data.  We  now  describe 
these  and  other  properties  which  it  may  be  desirable  for  / to  possess. 

2.1  Desirable  properties 

Monotonicity.  Here  we  require  / to  be  increasing  (respectively  decreasing)  if  (?/;)  is 
increasing  (respectively  decreasing).  More  generally  we  may  require  the  scheme  to  be 
‘co- monotone’,  i.e.  for  i = 0,  ...,N  — 1,  / is  increasing  (decreasing)  on  [xj,Xj+i]  if 
Vi  < Vi+i  (Vi  > Vi+ 1)-  Co-monotonicity  has  the  consequence  that  the  local  extrema  of 
/ occur  exactly  at  the  local  extrema  of  (yt).  Moreover  if  tji  = then  / is  constant 
on  [xj.Xi+i],  These  properties  may  be  too  restrictive  and  a weaker  alternative  is  what 
we  call  ‘local  monotonicity’:  for  i = 1, ....  A — 2,  / is  increasing  on  [xj, Xj+i]  if 
Pi-i  < Vi  < Vi+ 1 < yi+ 2 (and  similarly  for  decreasing).  Although  this  is  not  generally 
stated,  it  is  also  desirable  that  for  i = 0, . . . , A — 1,  / has  at  most  one  local  extremum 
on  (xi,Xi+i). 

Convexity.  Here  we  require  / to  be  convex  (concave)  if  the  piecewise  linear  interpolant  is 
convex  (concave).  More  generally  we  call  the  scheme  ‘co-convex’  if  for  i = 1, . . . , A - 2, 
/ is  convex  (concave)  on  [xj,Xi+ 1]  if  the  piecewise  linear  interpolant  is  convex  (concave) 
on  [xi-\,Xi+2]-  It  is  also  desirable  in  a co-convex  scheme  for  / to  have  at  most  one 
inflection  in  (xj,Xj+i),  0 < i < A - 1. 

Smoothness.  By  definition,  the  piecewise  linear  interpolant  is  shape-preserving,  and  so 
the  problem  is  trivial  unless  we  require  / to  have  greater  smoothness  than  continuity,  i.e. 
Ck  for  k > 1.  Since  all  the  schemes  use  piecewise  analytic  functions,  the  Ck  condition 
needs  to  be  checked  only  at  a finite  number  of  ‘knots’,  which  generally  include  the  data 
points.  We  remark  that  smoothness  and  shape-preservation  may  not  be  compatible;  e.g. 
if  for  i = 0, . . . ,4,  x*  = i — 2,  yi  = |xj|,  and  / is  convex  on  [xo,X4],  then  /(x)  = |x|, 
— 2 < x < 2,  and  so  is  not  C 1 at  0. 

Approximation  order.  It  is  generally  supposed  that  the  data  arise  as  values  of  some 
unknown  ‘smooth’  function  g,  i.e.  tji  — g(Xi),  i = 0, . . . , N.  Then  we  can  consider  how 
fast  the  interpolant  / converges  to  g as  we  increase  the  density  of  data  values  x,  in  the 
fixed  interval  [a, 6],  A scheme  has  approximation  order  0(hm)  if  ||/-</||  = 0(hm),  where 
h = max{xi+i  — x*  : i = 0, . . . , N — 1}  and  the  usual  norm  is  ||F||  = sup{|F(x)|  : a < 
x < b}. 

Locality.  In  a ‘global’  scheme,  the  value  /(x),  for  any  x,  generally  depends  on  all  the 
data.  In  contrast,  for  a ‘local’  scheme,  /(x)  depends  on  the  data  values  (x,,j/j)  only  for 
Xj  ‘near’  x.  There  may  be  advantages  in  local  schemes,  e.g.  when  data  are  modified  or 
inserted. 

Fairness.  It  is  often  desirable  that  the  curve  is  ‘fair’,  i.e.  pleasing  to  the  eye,  see  Section 

3. 
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Other  desirable  properties  are  invariance  under  scaling  or  reflection  in  x or  y.  and 
stability,  i.e.  small  changes  in  the  data  produce  small  changes  in  /.  There  may  also  be 
other  constraints  on  /,  e.g.  / > 0 when  j/i  > 0,  i = 0, . . . , N. 

2.2  Tension  methods 

Many  tension  methods  are  a modification  of  cubic  spline  interpolation,  which  we  now 
describe.  Given  data  (2.1),  there  is  a unique  function  / satisfying  (2.2),  where  / is  C2, 
is  a cubic  polynomial  on  [xi,Xi+ 1],  i = 0, — 1,  and  satisfies  suitable  boundary 
conditions  at  xo  and  xn . The  function  / minimises  f*"{g  )2  over  a suitable  class  of 
functions  and  this  energy  minimisation  property  is  generally  considered  to  give  a fair 
curve.  Determining  / requires  solving  a global,  strictly  diagonally  dominant  tridiagonal 
system  of  linear  equations. 

Since  cubic  spline  interpolation  is  not  shape-preserving,  in  1966  Schweikert  [67]  mod- 
ified the  scheme  by  replacing  cubic  polynomials  on  each  interval  [xt,Xi+ 1]  by  solutions 
of 

/(4)  - A if"  = 0, 

where  Aj  > 0.  When  Aj  = 0,  / will  reduce  to  a cubic,  while  as  A;  — > oo,  / approaches 
a linear  polynomial.  Thus  A*  acts  as  a tension  parameter  and  by  making  appropriate 
choices  of  Aj  large  enough  the  function  will  preserve  monotonicity  and/or  convexity 
globally  or  locally. 

Many  papers  have  been  written  on  Schweikert’s  tension  splines  giving,  for  example, 
ways  of  choosing  the  values  of  the  tension  parameters,  e.g.  [68,57,46,60].  However  the  fact 
that  the  method  uses  exponential  functions  can  be  seen  as  a drawback.  An  alternative  was 
introduced  by  Nielson  in  1974  [55]  by  adjusting  the  minimisation  property  of  cubic  splines 
to  a minimisation  problem  involving  also  the  first,  derivative.  The  resulting  function, 
called  a r'-spline,  is  also  cubic  on  each  interval  \xl,xl+\)  but  only  C1.  However  the  form 
of  the  C1  continuity  gives  extra  ‘smoothness’  for  parametrically  defined  curves  and  so 
we  discuss  z/-splines  further  in  Section  3.  By  generalising  the  minimisation  problem  still 
further  one  can  gain  a C1  piecewise  cubic  interpolant  with  further  parameters  for  gaining 
shape  properties  [22], 

The  idea  of  using  rational  functions  in  tension  methods  was  introduced  by  Spath  [69], 
also  in  1974,  and  put  in  a general  setting  of  tension  methods  in  [57].  Prom  1982-1988, 
Gregory  and/or  Delbourgo  produced  a series  of  algorithms  using  rational  functions,  e.g. 
[19,36,20,21,18].  We  illustrate  the  ideas  with  an  algorithm  from  [37].  Here  / is  C2  and 
on  each  interval  [xi,Xi+\]  it  has  the  form,  for  some  a,  b,  c,  d, 

„ , . _ a + bt  + ct2  + dt3  _ x - Xj 
1 + Ajt(l  — t)  ’ Xi+i  — Xi 

For  Aj  > — 1,  i = 0, . . . ,N-l,  f can  be  determined  as  the  solution  of  a strictly  diagonally 
dominant  tridiagonal  linear  system  (and  hence  the  scheme  is  global).  When  all  Aj  = 0, 
/ reduces  to  the  usual  cubic  spline  interpolant,  while  as  A;  — + oo,  / converges  uniformly 
to  the  linear  interpolant  on  [sj,a;j+i].  In  general  the  approximation  order  is  0(h2)  for 
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data  from  a C4  function.  In  the  special  case  of  monotone  data,  choosing 

\i  = Vi  + (f' fa)  + f (xi+ i))^±i — — , Hi>-  3,  i = 0, . . . ,1V  — 1, 

2/i+i  - 3/; 

ensures  that  / is  correspondingly  monotone,  and  for  the  choice  p,  = —2,  / reduces  to  a 
rational  quadratic  which  gives  optimal  approximation  order  0(h4).  Similarly  for  convex 
data,  / is  also  convex  provided  that  each  A,  satisfies  an  inequality  involving  / (xi), 
f (£j+i),  and  choosing  A;  appropriately  (which  requires  solving  a non-linear  equation) 
further  ensures  approximation  order  0(h 4). 

There  are  some  more  recent  methods  involving  rationals,  e.g.  [58]. 

The  idea  of  using  variable  degree  to  preserve  shape  was  introduced  by  McAllister, 
Passow  and  Roulier  in  1977  [47,56].  They  produce  monotone,  convex  schemes  of  arbit- 
rarily high  smoothness  by  constructing  a shape-preserving  piecewise  linear  interpolant  l 
with  one  knot  between  any  two  data  points  (and  no  knots  at  the  data  points)  and  then 
defining  the  final  interpolant  on  each  interval  [xi,xi+ 1]  as  the  Bernstein  polynomial  of  l 
of  some  degree  m,.  The  idea  was  extended  from  1986  by  Costantini  [8-10].  For  k > 1, 
m,'  > 2k  + 1,  i = 0, . . . , N — 1,  he  constructs  a shape-preserving  piecewise  linear  inter- 
polant l with  knots  at  Xi  + k(xi+i  -Xi)/m,i  and  a;,+i  - k(xi+ \ -Xj)/mj,  i = 0, . . . ,N—  1. 
The  final  interpolant  / coincides  on  each  interval  [rr,,a;,+i]  with  the  Bernstein  polyno- 
mial of  l of  degree  m,  and  is  hence  Ck  (with  f!ji{xi)  — 0,  j = 2, . . . , k).  In  [10]  there  is 
a co-monotone,  co-convex  scheme  in  which  the  degrees  m,  can  either  be  chosen  a priori 
or  computed  automatically  according  to  the  data. 

The  above  schemes  using  variable  degree  are  not  strictly  tension  schemes  in  our  sense 
but  in  1990,  Kaklis  and  Pandelis  [40]  introduced  a tension  method  by  using  the  above 
form  for  k — 1,  i.e.  on  each  interval  [r,;,  r,+i]  it  has  the  form: 

f{t)  = f(xi)(  1 - t)  + f(xi+i)t  + Cit(  1 - t)mi  + ditmi(  1 - f),  t = — — — . 

Xi+i  - Xi 

Here  rrij  > 2 is  an  integer  and  for  each  choice  of  mo, . . . ,mjv_ i,  the  numbers  c,,  d.t  are 
chosen  so  that  / is  C2,  which  requires  the  solution  of  a strictly  diagonally  dominant 
tridiagonal  linear  system.  When  all  m,  = 2,  this  reduces  to  the  usual  cubic  spline  inter- 
polant, while  as  rrii  — > oo,  / converges  uniformly  to  the  linear  interpolant  on  [xi,  r,+i] 
with  order  Ofm”1)  (or  0(m~ 2)  if  m,_i,  m,;+1  remain  bounded).  For  further  discussion 
of  variable  degree  shape-preserving  functional  interpolation,  see  [11]. 

Our  final  type  of  tension  method  was  introduced  by  Manni  [50]  in  1996.  The  general 
idea  is  to  define  / on  [xi,Xi+i]  as 

f(x)  =Pi{q~1{x)), 

where  p*,  qt  are  cubic  polynomials  on  [xi,Xi+i]  and  Qi  is  strictly  increasing  from  [xi.xt+i\ 
onto  itself,  so  that  the  inverse  q~l  is  well-defined  on  [x,i, x.;+i].  For  / (xi)  = di,  i = 
0, . . . , N,  we  require 

Pi  (t  ) Ajdj,  qi{xi)  — Xi,  Pi{xi+ 1)  pidi+i,  q^  (x^_j_] ) = 

for  parameters  A,  > 0,  p,  > 0.  For  A,  = = 1,  we  have  qdx)  = x and  / reduces  to  a 

cubic  on  [xi,Xi+i],  while  for  Aj  = pi  = 0,  / becomes  linear  on  [xi,Xi+i]. 
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In  [50],  the  values  do, dm  are  assumed  known  (or  estimated  from  the  data  values) 
and  the  scheme  is  local  C 1 , gives  necessary  and  sufficient  conditions  for  the  values  of  the 
parameters  A,,  /i,  forco-monotonicity,  and  has  approximation  order  0(h2)  when  g is  C2 
and  generally  0(h4)  when  g is  C4. 

Manni  and  co-workers  have  written  a series  of  papers  using  the  same  idea,  [51,53,54]. 
For  example  in  [45],  the  values  d,:  are  not  assumed  given  but  are  chosen  to  ensure  that 
the  function  is  C2,  thus  providing  a locally  monotone,  co-convex  global  scheme  which 
generalises  usual  cubic  spline  interpolation;  while  in  [52]  two  further  knots  are  inserted 
in  each  interval  [x,,Xi+ 1]  to  produce  a C2,  locally  monotone,  co-convex  local  scheme 
which  interpolates  values  of  /^(x?),  j = 1, 2,  i = 0, . . . , N. 

2.3  Direct  methods 

In  1967,  Young  [71]  considered  shape-preserving  interpolation  by  polynomials  and  a num- 
ber of  papers  have  appeared  since  on  this  topic,  e.g.  [59]  gives  a constructive  proof  of 
the  existence  of  a co-monotone  interpolant  with  an  upper  bound  on  the  degree  required. 
However  for  a practical  algorithm,  using  a piecewise  polynomial  offers  much  more  flexib- 
ility than  a single  polynomial.  Numerous  papers  have  been  written  using  such  polynomial 
splines  and  we  mention  briefly  only  a few. 

By  inserting  extra  knots  between  data  points,  a convexity  preserving  scheme  with 
C2  cubics  was  given  by  de  Boor  [4,  p.303],  and  co-monotone,  co-convex  schemes  with 
C 1 quadratics  in  [48,49,66].  C1  cubic  splines  with  knots  at  the  data  points  are  used  for 
co-monotonicity  in  [25,5,24,70],  (the  last  of  these  using  a variational  approach),  and  for 
both  co-monotonicity  and  co-convexity  in  [16,17].  We  also  recall  the  methods  using  spline 
functions  of  variable  degree  with  knots  between  the  data  points  to  obtain  interpolants 
with  arbitrarily  high  smoothness  which  were  discussed  under  tension  methods. 

Finally  we  note  that  following  the  paper  [62]  which  was  as  early  as  1973,  Schaback 
[63]  gives  a C2  co-monotone,  co-convex  scheme  which  uses  a cubic  polynomial  on  any 
interval  [xi,Xi+ 1]  where  an  inflection  is  needed,  and  on  other  intervals  employs  a rational 
function  of  form  quadratic/linear. 

3 Planar  curves 

Given  data 

Ii€R2,  i = 0,...,N, 

we  consider  a curve  r : [a,  5]  — > R2  satisfying 

r(ti)  = Ii,  i = 0, . . . , N,  (3.1) 

for  values  a = to  < h < — b.  For  a closed  curve  the  situation  is  extended 

periodically  so  that 

Ii+N  = Io,  ti+N=ti , iez,  r(t  + b — a)  = r(t),  t E R. 

3.1  Desirable  properties 

Shape.  For  this  case  it  is  not  usually  relevant  to  consider  preservation  of  monotonicity. 
We  say  a scheme  is  ‘co-convex’  if  the  curve  r has  the  minimum  number  of  inflections 
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consistent  with  the  data.  In  practice,  schemes  satisfy  the  somewhat  stronger  condition 
that  for  any  0 < i < j — 2 < N - 2,  r is  positively  (negatively)  locally  convex  on 
[ti+i,tj-i]  if  the  polygonal  arc  joining  Ij  is  positively  (negatively)  locally  convex. 

For  more  details  on  this  and  other  desirable  properties,  see  [29]. 

Smoothness.  We  shall  call  the  interpolating  curve  Ck  for  k > 0 if  the  function  r is 
Ck.  A G°  curve  r we  shall  call  G1  if  the  unit  tangent  vector  is  continuous,  and  G2  if,  in 
addition,  the  curvature  is  continuous.  A Ck  curve  r is  Gk,  k = 1,2,  provided  that  the 
parameterisation  is  regular,  i.e.  r ( t ) ^ (0,0),  which  is  generally  desirable.  It  is  usually 
sufficient  to  have  Gk,  rather  than  Gfc,  continuity  if  only  the  appearance  of  the  curve  is 
important  and  the  choice  of  parameter  t is  not  significant. 

Fairness.  Planar  curves  often  arise  in  computer-aided  design  where  it  may  be  par- 
ticularly important  that  the  curve  is  pleasing  to  the  eye.  Though  this  is  subjective, 
various  criteria  have  been  suggested  to  be  relevant,  such  as  magnitude,  rate  of  change 
or  monotonicity  of  the  curvature.  Some  schemes  include  ‘shape  parameters’  which  can 
be  manipulated  by  the  designer  to  modify  the  shape  of  the  curve. 


Approximation  order  is  not  important  in  the  context  of  design  when  the  data  are 
not  considered  to  be  taken  from  some  unknown  curve.  Approximation  order  is  related 
to  reproduction  of  polynomial  curves,  and  a related  property  for  planar  curves  is  re- 
production of  arcs  of  circles  (or  more  generally  conics);  this  cannot  be  done  exactly  by 
polynomials  but  it  can  be  achieved  by  using  rationals. 

Locality  and  other  desirable  properties  are  similar  to  the  functional  case  as  described 
in  Section  2.1,  though  it  is  generally  more  appropriate  that  the  invariance  is  under  a 
rotation  and  the  same  scaling  in  both  x and  y. 


3.2  Tension  methods 


In  Section  2.2  we  mentioned  Nielson’s  ^-splines  [55].  Applying  this  scheme  for  both 
components  of  r gives  a function  r which  is  cubic  on  each  interval  [ti,  tj+i],  is  G1  and 
satisfies 


r"  (ft)  = r"  (tj  ) + Vir\ti),  i = 1, . . . ,N  - 1, 


where  Ui  > 0.  This  condition  is  sufficient  for  G2  continuity  of  r (assuming  regular  para- 
meterisation). When  all  Vi  = 0,  r will  reduce  to  the  usual  G2  cubic  spline  interpolant.  As 
Ui  — * oo,  the  curve  is  ‘pulled  tight’  at  fy  and  as  Ui,  Ui+i  — > oo,  it  approaches  the  linear 
interpolant  on  [ti,tj+i]. 

The  scheme  in  [37]  by  Gregory  which  was  mentioned  in  Section  2.2  was  adapted  to 
the  planar  case  in  [38],  Other  schemes  using  rationals  were  proposed  by  Clements  in 
[6,7],  where  r is  a G2  curve  which  on  each  interval  [U,U+ 1]  has  the  form,  for  some  a,  b, 
c,  d £ R2, 


. . a(l  — s)3  , . 

r(f)  = — — + 6(1  - s)  + cs  + 

WiS  + 1 


ds3 

Wi(\  - s)  + 1! 


t-tj 
t'i+1  ti 


where  Wi  > 0 are  the  tension  parameters. 
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The  variable  degree  tension  method  of  [40],  also  mentioned  in  Section  2.2,  was  adapted 
to  the  planar  case  in  [41],  and  extended  in  [27]  to  allow  the  designer  to  obtain  a ‘fair’ 
curve  by  minimising  the  number  of  changes  in  the  monotonicity  of  the  curvature. 

3.3  Direct  methods 

The  papers  [34,35,28,23]  give  local,  G 2 eo-convex  schemes,  e.g.  in  [28],  a rational  cu- 
bic/cubic is  used  on  each  interval  [tj,ti+i]  and  the  tangent  vectors  and  curvatures  are 
stipulated  by  the  algorithm  to  ensure  that  the  convexity  conditions  are  satisfied  and 
circular  arcs  are  reproduced,  with  the  possibility  of  modifying  the  tangent  vectors  and 
curvatures  further  as  shape  parameters. 

Following  an  earlier  scheme  in  [64],  Schaback  in  [65]  gives  a global  G 2 co-convex 
scheme  which  uses  a cubic  polynomial  on  any  interval  [U,U+ j]  where  an  inflection  is 
needed,  and  on  other  intervals  employs  quadratic  polynomials. 

Sapidis  and  Kaklis  [61]  give  a G2  co-convex  scheme  by  interpolating  by  a piecewise 
quintic  curve  tangent  directions  and  curvatures  gained  by  their  tension  method  [41]. 

In  [1]  a local,  co-convex  G2  scheme  is  given  which  uses  polynomials  of  degree  six 
and  which  attempts  to  obtain  a fair  curve  by  imposing  conditions  on  the  curvature  to 
minimise  measures  of  fairness.  Finally  we  note  that  in  [12]  Costantini  gives  an  abstract 
theory  and  general  purpose  code. 

4 Space  curves 

Given  data 

ReR3,  i = 0,...,JV, 

we  consider  a curve  r : [ a,b } — > R 3 satisfying  condition  (3.1)  as  before. 

4.1  Desirable  properties 

What  is  meant  by  ‘shape-preserving’  is  not  so  clear  for  space  curves  as  for  the  planar 
case.  Criteria  were  introduced  by  Kaklis  and  Karavelas  [39]  and  extended  by  Ong  and 
the  author  in  [31].  We  shall  sketch  these  below.  They  are  discussed  in  further  detail  in 
[30],  where  some  further  extensions  are  suggested.  We  write,  for  appropriate  indices  i: 

Li  — /j,  Aj  = det [Tj— i , Z/j, Ti-i-i],  Nj  = Li—\  x L(. 

Torsion.  We  ensure  that  the  curve  is  ‘twisting’  in  the  same  manner  as  the  piecewise 
linear  interpolant  by  requiring  that  if  A,-  ^ 0,  then  the  torsion  of  r has  the  same  sign  as 
A,  on  ( [tii 

Convexity.  Let 

K(t)  = r (t)  xr  (f),  a <t  <b. 

We  require  that  for  1 < i < N — 1,  K(ti).Ni  > 0,  which  means  that  the  projection  of 
the  curve  r onto  the  plane  of  R-i,  L,  R+i,  has  the  same  sign  of  local  convexity  at  as 
the  polygonal  arc  Ii-iRIi+i.  Moreover  if  A^.W+i  > 0,  we  require 

K(t).Nj  > 0,  j = i,i  + 1,  ti<t<ti+ 1, 
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which  implies  that  the  curve  r has  the  same  sign  of  local  convexity  on  when 

projected  in  any  direction  AATj  + (1  — X)Ni+i  for  0 < A < 1.  Finally  we  require  that 
if  Ni.Ni+i  < 0,  then  for  j = i,  i + 1,  K{t).Nj  has  exactly  one  sign  change  in 
which  imples  that  each  of  the  above  projections  of  r have  just  one  inflection. 

Smoothness.  This  is  as  for  planar  curves,  except  that  we  call  the  curve  G3  if  it  is  G2 
and,  in  addition,  the  torsion  is  continuous.  Other  desirable  properties  are  similar  to  the 
planar  case. 

4.2  Tension  methods 

Although  interpolation  by  space  curves  with  a special  shape  is  considered  in  [44],  the 
first  specific  shape-preserving  interpolation  scheme  by  space  curves  was  due  to  Kaklis 
and  Karavelas  [39],  who  adapted  the  variable  degree  tension  method  of  [40]  to  give  a 
G2  method  which  was  also  G3,  but  at  the  expense  of  zero  torsion  at  the  data  points.  In 
[42]  the  same  authors  adapted  Nielson’s  ^-splines  to  the  three  dimensional  case  to  give  a 
curve  which  is  C1  and  G2.  The  paper  [14]  also  uses  variable  degree  for  tension  parameters 
but  gives  a G3  scheme  in  which  the  limiting  curve  as  the  tension  goes  to  infinity  is  not 
the  piecewise  linear  interpolant  but  the  shape-preserving  interpolant  given  by  either  of 
the  above  two  schemes.  In  [15]  a G3  scheme  is  also  given  but  here  the  components  of  r 
on  each  interval  [f,,fi+i]  lie  in  the  linear  span  of  the  functions 

(1  - u),u,  (1  - u)*\  (1  - u)mi-\  (1  - u = 1 ~ ■ *■■■  . 

U+l  ~ H 

When  mj  = mi+i  = 5,  this  reduces  to  a quintic  polynomial.  As  m*,  mi+\  — > oo,  it  tends 
to  a linear  polynomial  and  then  the  curve  r approaches  the  piecewise  linear  interpolant 
on  [ti,ti+i\. 

The  paper  [26]  also  uses  variable  degree  splines  with  degree  on  each  interval  at  least 
five,  and  the  curve  r also  converges  to  the  piecewise  linear  interpolant  as  the  degrees  go 
to  infinity.  However  here  the  curve  is  C4,  which  the  authors  feel  may  give  extra  fairness 
to  the  curve  due,  for  example,  to  lowering  the  maximum  absolute  value  of  the  curvature. 
Variable  degree  polynomial  splines  are  also  used  in  [13]. 

4.3  Direct  methods 

Following  an  earlier  scheme  in  [31],  Ong  and  the  author  gave  a local  G2  scheme  in  [32] 
which  employed  a rational  cubic/cubic  between  data  points,  extending  the  ideas  of  the 
planar  scheme  in  [28].  This  was  further  extended  to  a local  G3  scheme  using  a rational 
quartic/quartic  in  [43].  In  [33],  the  degrees  of  freedom  inherent  in  the  scheme  in  [32] 
were  used  to  optimise  a fairness  measure.  Finally  we  mention  the  papers  [2,3]  which  give 
local  G3  schemes  using  a piecewise  polynomial  of  degree  six,  also  allowing  optimisation 
of  a fairness  measure. 

It  will  be  noted  that  many  of  the  above  papers  are  extremely  recent  and  it  is  hoped 
that  the  unavoidable  lack  of  detail  here  will  serve  to  tantalise  readers  to  discover  for 
themselves  more  of  this  rapidly  developing  field. 


T.  N.  T.  Goodm,an 
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